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^ ■ Abstract 

^ . This paper is a theoretical work in support of a newly proposed experiment (R. Stutz and E. 



Cornell, Bull. Am. Soc. Phys. 89, 76 2004) that promises greater sensitivity to measurements 
of the electron's electric dipole moment (EDM) based on the trapping of molecular ions. Such 



K^ ■ an experiment requires the choice of a suitable molecule that is both experimentally feasible and 

in ; 

^^ . possesses an expectation of a reasonable EDM signal. We find that the molecular ions PtH^, 

^T \ HfH+, and HfF+ are suitable candidates in their low-lying ^A states. In particular, we anticipate 



that the effective electric fields generated inside these molecules are approximately of 73 GV/cm, 



Y^ • -17 GV/cm, and -18 GV/cm respectively. As a byproduct of this discussion, we also explain how 

^ \ to make estimates of the size of the effective electric field acting in a molecule, using commercially 



available nonrelativistic molecular structure software. 
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I. INTRODUCTION 

The Standard Model of elementary particle physics admits that the electron may have 
an electric dipole moment (EDM), but a darn small one; the Standard Model prediction is 
on the order of d^, = 10~^^ e-cm, many orders of magnitude below the current experimental 
limit. Nevertheless, in models beyond the Standard Model, notably super-symmetry, a larger 
EDM is expected, sometimes tantalizingly close to the current experimental upper bounds 
on de |2j. To observe dg would imply a different form of charge-conjugation/parity (CP) 
violation than the familiar form in the K- and B-meson systems. The observation of de 
requires extremely large electric fields to observe the energy shift due to different alignments 
of de along the quantization axis. 

The current upper limit on the electron EDM, 1.6 x 10^^^ e-cm, comes from measurements 
in atomic thallium [jj. In this experiment, an effective electric field is generated inside the 
atom. The magnitude of this electric field (on the order of 70 MV/cm) is far larger than 
the field that could be applied directly in the laboratory. This circumstance is what makes 
an atomic EDM experiment desirable. 

Taking this idea even further, Sandars noted that the effective electric field within a 

n 

molecule can be even larger than that within an atom j3|. This insight has prompted a 
series of proposed and actual experiments, aimed at using molecules as high-electric-field 
laboratories. These experiments have _emgloyed either ^molecular beams or traps for neutral 

30 14 |5| and YbF |8l |9||. Several others have also been 



molecules, including molecules PI 

p..oposed a. viable candidates BflaflQQ. 



Recently, a new experiment has been proposed that may increase the experimental sen- 
sitivity to an electron EDM by orders of magnitude [ij, |l5[ . The key to this experiment is 
using molecular ions in lieu of neutral molecules. Because ions are far easier to trap for long 
times than neutrals, the proposed experiment wins in terms of long coherence times. Never- 
theless, such an experiment poses a number of challenges, both in experimental technology 
and in the choice of a reasonable molecule. 

The purpose of this paper is to identify diatomic molecular ions that possess properties 
that are desirable for this experiment. For instance, the molecule should have a relatively 
small degeneracy of quantum mechanical states near its ground state. This feature, in a 
cold sample at thermal equilibrium, would ensure that most or all of the ions are already in 



or near a state that is useful for the EDM measurement. Because the number of ions in a 
trap is fairly small, this requirement is necessary for achieving a good signal-to-noise ratio. 
To meet this requirement we restrict our attention to diatomic hydrides, where the light 
hydrogen atom will contribute to a relatively large rotational constant. Similarly, we are 
interested in heavy atoms with spinless isotopes in order to ensure a favorable population 
distribution. This requirement is an experimental one. On the theory side, isotopes with 
spin are useful in characterizing the accuracy of the calculations and for making empirical 
estimates |^, |l3| . A second requirement is that the electric field used to polarize the molecule 
must be relatively small, on the order of 100 V/cm or less. This is because in the proposed 
experiment, the polarizing field must be rotating at MHz frequencies in order to act on 
the molecules without accelerating them out of the trap. Lastly, we restrict ourselves to 
molecules that have at least one unpaired spin (paramagnetic, i.e. molecules with S ^ 0), 
so that the net EDM of the electrons does not vanish 

Because we confine our attention to diatomic molecules, we must meet two additional 
requirements: (1) A large dipole moment, on the order of 1 Debye, which is not hard to 
achieve if the two atoms are quite different, and (2) A small energy splitting between even- 
and odd-parity branches of the ground state. In other words, a small omega doublet is 
desirable. 

In general, the size of the omega doublet is smaller for larger quantum number fl, which 
represents the total (orbital plus spin) angular momentum of the electron about the inter- 
molecular axis. For this reason we prefer diatomic molecules with high angular momentum. 
In the following we consider species with ^A molecular symmetry, which can support values 
of Q as large as f2 = 3. This choice is in sharp contrast with previous EDM candidate 
molecules, most of which were based on ^S molecular symmetries. More recently, molecular 
ions of ^E y, |5[ or ^11 [ij, |l3[ symmetry have been considered as well. 

Finally, the essential requirement for the molecule is that it must possess a large internal 
effective electric field. It is this field that acts on the electron's EDM to make a measur- 
able signal. A quick review of the relevant physics reveals why a large internal field is so 
important. If an electric field were applied to an atom, and if relativity did not matter, 
then a given electron would feel no net electric field. In equilibrium, the applied external 
field would be exactly balanced by the field due to all the other charges in the atom, or 
else the electron would be pulled to a different equilibrium. This is the content of Schiff 's 



theorem [16| . However, the fact that the electron can move at relativistic speeds in an atom 
provides a loophole around this theorem. In the relativistic case, the electron's equilibrium 
is governed by the balance of electric fields, plus the motional magnetic field generated when 
it swings by the nucleus. Thus the net electric field, by itself, can be non vanishing. By 
applying an external electric field, the non zero field within the molecule can be increased 
by a relativistic enhancement factor manifested through the mixing of s and p atomic states. 
The resulting energy shift scales with the third power of the nuclear charge, and therefore 
strongly favors heavy atoms as EDM candidates. 

Qualitatively, molecules can be even more effective at mixing s and p states of atoms than 
externally applied fields. Indeed, in molecular orbitals it is possible for s and p levels to 
perfectly hybridize, as in the conjugated bonds of carbon-bearing molecules. The relativistic 
enhancement factor can be estimated from its value for the heavy atom, therefore estimating 
the effective internal electric field becomes a question for nonrelativistic molecular structure 
theory. The mixing of s and p atomic orbitals is useful only when it produces a molecular 
orbital of a symmetry. An electron in such a molecular orbital has no orbital angular 
momentum about the molecular axis, and can thus penetrate close to the heavy atom's 
nucleus, where relativistic effects are greatest. This behavior explains why molecules with 
total symmetry S have been preferred in the past. 

As noted above, however, the experiment under consideration is likely constrained to 
molecules with nonzero angular momentum about the axis. To accomplish these seemingly 
contradictory goals, we seek molecules having two valence electrons. One, of a molecular 
symmetry, will be responsible for the effective electric field and consequent EDM signal. The 
second, of 6 symmetry, will allow for small omega doubling. We are thus led to contemplate 
molecules of {a6)^AQ symmetry overall. Further, since 6 molecular orbitals are likely to 
arise from d-tjpe atomic orbitals, we seek diatomic molecules in which the heavy atom is a 
transition metal. Previous theoretical work on heavy transition metal hydrides |l2| identifies 
a pair of likely molecules with low- lying '^A states, namely, HfH"'" and PtII+. 

In the following sections we will explore the properties of these and other molecules in 
terms of their usefulness to an electron EDM experiment. In Sec.|n]we revisit the theoretical 
methods for estimating the EDM signal in terms of nonrelativistic molecular wave functions 
and the way in which such wave functions can be constructed. In Sec. lIIII we present concrete 
examples of calculations for various diatomic molecules. These calculations verify that the 



perturbative methods we employ can semi quantitatively reproduce the effective electric 
fields predicted previously by more substantial calculations. 

II. BACKGROUND: OBSERVABLE CHARACTERISTICS OF THE ELECTRON 
EDM 

A. The effective electric dipole moment of an atom 

In EDM circles, it is customary to speak of an "effective electric field" experienced by the 
electron inside an atom or molecule. What this means in practice is that the energy shift 
due to an electron EDM of magnitude de is simply 

A£; = -4Feff. (1) 

We use the letter F to denote electric field strengths, to distinguish them from energies, 
denoted E. To evaluate the energy shift, one writes the perturbation Hamiltonian in Dirac 
notation ([211, p. 254 ff.) is: 

\^0 24a -Fy 
This 2x2 matrix acts on the state vector ip = {ipjji V'l)? where ipu and ipL are the "upper" 
and "lower" components of the electron's Dirac wave function. Each oiipu and ipL-, in turn, 
stands for a two-component spinor. 

In Eq.(j2I), 0" is a Pauli matrix designating the electron's spin, so the combination df,(T 
describes the EDM of the electron. (The direction of the electron's EDM is presumed to 
be collinear with its magnetic moment). This moment interacts with the instantaneous 
electric field F experienced by the electron. To a good approximation, this field is given 
by the Coulomb field of the heavy nucleus with charge Ze, F = Zefjr^. Here r is the 
electron-heavy atom distance. Here and throughout this paper we use atomic units, setting 
^ = me = 1. 

Next we denote the lower component of the electron's ground state wave function to be 
1'?/'°). Then the perturbative energy shift due to H^ would be ('?/'° |2(ieO"-F|'?/'£,)- However, this 
expression vanishes for a ground-state wave function of definite parity, owing to the odd- 
parity of F. Thus, an externally applied electric field is required to mix states of different 



parities. In the important example of an atom with an s-electron (e.g., Cs), an apphed 
electric field would create a linear combination of s and p states, 

K) = es\Ms)) + e,\Mp))- (3) 

Higher angular momentum states are also possible in principle but do not make a substantial 
contribution at laboratory fields. In typical atomic EDM experiments, the field is sufficiently 
weak to invoke the perturbative approximations: 

e, ^ 1 (4) 

(s|er-Fappiicdb) /gv 



e. 



r-^ 



'" E{s) - E{p) 



Here er stands for the relevant dipole operator, and the denominator denotes the energy 
difference between the s and p states. 

With these approximations, the expression for the energy shift becomes 

m2d,a ■ F\i;l) = 2e.e, (^) r,ei. (6) 

The factor Frei incorporates the relativistic effects and can be estimated analytically using 
relativistic Coulomb wave functions 

where a is the fine structure constant and 7 = a/(j + 1/2)^ — (Za)^. The quantities Ug 
and Up are the effective quantum numbers for the s and p states of the heavy atom, defined 
as the the principal quantum number shifted by a quantum defect. All values of Ug and 
Up are obtained from Ref. 43|. The constant Zi is the effective nuclear charge seen by the 
valence electron; for a neutral atom Z^ = 1. An additional energy shift would arise from 
perturbations due to the P3/2 excited state. However, the matrix element connecting the 
..,, state to the p„. state ,s van,shi„g.y small, smce the p,,_ radial funcfo,. ha. „egl,g,ble 
amplitude near the nucleus |21| . 

The energy shift of an atom in an electric field, due to the EDM of the electron thus 
requires mixing the S1/2 ground state and the pi/2 excited state. The resulting energy shift 



IS, m atomic units. 



AE 



4(si/2kbi/2) ^ Ze 

-J- rel 



E{si/2) - E{pi/2) '^ al 



2" 

C^e-^applied- (8) 



This expression can be viewed in two ways. One way is to envision that the atom gains 
an EDM that is bigger than the electron's EDM by the factor in square brackets in Eq.(jH)). 
Alternatively, the existing electronic EDM may be seen as experiencing the applied field 
-^applied) enhanced by this same factor. 

B. The effective electric field inside a molecule 

In diatomic molecules, the same kind of s-p mixing may occur, perhaps much more 



substantially than in an atom |22l. |23|. If so, the basic expression for the measured energy 
shift is still given by Eq.®, but the coefficients may have quite different values. To estimate 
these coefficients, we expand the molecular orbital wave function of the contributing electron 
in terms of atomic orbitals: 

|mol) = e^|s) + ep\p) + ^ Cotherl other). (9) 

other 

Here s and p still refer to the relevant atomic orbitals on the heavy atom, from which the 
EDM shift will arise. But now there are other orbitals that may participate. For example, 
the electron wave function may have considerable contributions from d or f orbitals on the 
heavy atom or from atomic orbitals residing on the other atom. Thus, in general, it may not 
be a good approximation that ef + ep = 1. Nevertheless, one may seek molecules for which 
|e^| ^ [cpl ^ l/v2, which would be the optimal mixing for our purposes. To determine these 
quantities for a given molecule is a problem in nonrelativistic molecular orbital theory. 

Once the amplitudes e^ and e^ are known, we can use Eq.ijHj) to find the energy shift, with 
the following caveat. The orbital \p) here stands for the nonrelativistic p wave function of 
the atom. It is therefore a linear combination of the relativistic atomic states, denoted by 
total spin j: 

b) = -^bi/2) + y 3^3/2). (10) 

The factor 2a emphasizes that the pi/2 contributes with opposite signs for the two spin 
projections on the molecular axis, a = ±1/2. As noted above, the P3/2 orbital does not 
contribute to the electron EDM shift; hence the relevant matrix element reads 

{s\2d,a ■ F\p) = --^(si/2|24a ■ FIP1/2). (11) 



Thus, following Flambaum 2M ^^^ Khriplovicli and Lamoreaux 23[, we find the electron 



EDM energy shift for a diatomic molecule (in atomic units) 

AE = 



4a Ze' 

. — (^s^^p^ rel 



^. (12) 

eao 



V3 ao 

In strong contrast to Eq.®, the molecular expression Eq.(fT^ does not depend explicitly 
on the applied electric field FappUod- The effective electric field is simply the coefficient in 
front of de- Applying an external field in the laboratory is essential, however. In zero electric 
field, the molecule's energy eigenstates are also eigenstates of parity ^3(1 |: 

\\Q\,±) = ^m±\-Q)), (13) 

but separated by an energy Edonh referred to as the f2-doubling energy. By convention, the 
lower (upper) sign of this doublet is denoted by the letter e (/). In such a state, the electronic 
spin projection a will lie along the molecular dipole moment as often as against it, canceling 
the energy shift. However, in a modest electric field FappHed, the signed values of ±Q are 
again good quantum numbers, and the expression Eq. (fT^ is relevant. To overcome -Edoub, 
the field must be large compared to a "critical" field, at which the Stark energy overcomes 
the energy difference between the two parity eigenstates: 

Fcrit = ^, (14) 

where (Im is the permanent electric dipole moment of the molecule. For the proposed 
molecular ion experiment, it is desirable to make Fait as small as possible, preferably below 
100 V/cm []J. 

III. COMPUTATIONAL METHOD 

There remains the question of evaluating the constants e^ and ep for a given molecule of 
interest. Fortunately, it is common practice in molecular structure theory to cast molecular 
orbitals as linear combinations of atomic orbitals (LCAO). The atomic orbitals, in turn, are 
usefully expressed in terms of a Gaussian basis set that is carefully constructed to reproduce 
the energy levels of the individual atoms. After decades of development, these methods 



are now robust and commercially available in software packages such as Gaussian 



MOLPRO 



2J and 



We use MOLPRO in the calculations performed here. 
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A. Molecular electronic structure 

We assume that the molecules are well-approximated by Hund's angular momentum cou- 
pling case (a). This approximation is justified by the much greater strength of the spin-orbit 
interaction and the exchange splitting (both on the order of lO'^cm"^) as compared to rota- 
tional energies (on the order of 10 cm~^). In this case the relevant quantum numbers of each 
electron are its projection of orbital (A) and spin (cr) angular momentum on the molecular 
axis. (In fact, strong spin-orbit mixing will impact these quantum numbers slightly, pushing 
the molecules toward Hund's case (c). We will deal with this later). 

To make a concrete calculation, each molecular orbital is expanded in a basis set of 
Gaussian functions, where a basis element is defined as 

\bj) = ^^i\9i)- (15) 

i 

Here each \gi) is a Gaussian function centered on either of the atoms in the molecule. Such 

a set of Gaussians is chosen to optimize the energies of each atom separately and is provided 

in one of a set of standard basis sets, referred to below. With this basis, it is now possible 

to define an orbital as 

\<Pk) = J2c,\b,), (16) 

j 
and a configuration is a product of occupations of orbitals given as 

\<^i)=Al[dk\<l)k). (17) 

k 

The dk can take the values 0,1 or 2. A stands for anti-symmetrization. 

The Hartree-Fock (HF) procedure minimizes the energy of the ground state configuration, 
hence only one configuration is optimized. This is done by minimizing the coefficients in the 
orbitals and configurations. For concreteness, the configuration is 

k j 

and the procedure finds the dk and Ckj that minimize the total electronic energy E. After 
this procedure, we perform a Mult i- Configuration Self Consistent Field (MCSCF) J2^ |22| 
calculation with a Complete Active Space (CAS) within MOLPRO. The MCSCF takes the 
HF E as an initial guess and finds appropriate linear combinations of the configurations 



given in Eq.()18p that minimize E over M configurations. The M configurations can contain 
many different symmetries and spin states. This wave function is 

I 

I k j 

and the procedure varies //, dik, and Ckj in order to minimize the average E. This is the 
time-consuming part of the calculation due to the non linear optimization routine. This 
step is then followed by the Multi-Reference Configuration Interaction (MRCI) 2^ |29| 
program. A MRCI calculation uses the MCSCF configurations and mixes in single and 
double excitations. The wave function is 

I^MRCI^^^^^I^MCSCF^^ (20) 

rn 

where this procedure finds the h^ that minimize the value of E for a given symmetry, subject 
to the constraint that the fi, dik and Ckj are all fixed. 

Therefore, deep inside the MCSCF wave function is the a molecular orbital that mixes the 
atomic s and p orbitals. We denote this orbital as |mol, a), which has an explicit expansion 
into the Gaussian basis functions according to Eq. (jl8|) . Analogously, the wave function for 
the heavy atom (e.g., Hf^ or Pt+) is also computed at the MCSCF and MRCI levels and 
contains atomic orbitals |atom, s) and |atom, p) expanded into the same Gaussian basis set. 
It is now a straightforward to estimate the amplitudes es,p = (mol|atom, {s,p)). 

There are two final, but sometimes decisive, corrections. First, the desired configuration 
does not constitute the entire wave function, but only a fraction of it, given by its amphtude, 
e.g. ho in the MRCI expansion in Eq. (j2(Jj) . For the molecules considered in the proposed 
experiment, however, we have found this factor to be within several percent of unity. Alter- 
natively, it does play a crucial role in the case of PbO where the factor ho is considerably 
less than unity for the (cra*)^S metastable state; the state that contributes to the EDM 
signal. 

More significantly, the molecular state identified in Eq. ()2()j) may be mixed with others 
via spin-orbit couplings. In the case that one other molecular symmetry is mixed in, we get 

1^) =cosx|^^+^AS,f])+sinx|^^'+^A'S',f]), (21) 
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in terms of a mixing angle x- There is no prime on the Vt since spin-orbit interactions preserve 
the value of Vt. The factor cosy can be quite significant. For example, in a ^11 state (as in 
PbF H a„d the Halide. HB.^ U a.d HI^ QB). the s.g.e va.e„ee ele*o„ would be ,n 

a 7r state and would not contribute at all to an effective electric field (cosx = 0). However, 
spin-orbit mixing with nearby ^S states can introduce a a electron, but its effective electric 
field is reduced by a factor of sin x- This mixing is responsible for the large effective field 
in PbF because it has a very large spin-orbit constant. Finally, we arrive at the expression 
used to estimate the effective electric field in our approximation: 



'cff 



A '7 

-^/iocosx(niol, 0" I atom, s)(mol, cr|atom, p)rrei— 3- 
V 3 o,Q 



(22) 



de 

in terms of analytical expressions and concrete features of MOLPRO output. 

B. Potential curves 

In the available literature, there is no consensus as to the true ground state of PtH"*" 
|l7L Il8l . Il9| : only one paper attempts to work out the ground state of HfH+ jl2|. The 



debate arises in PtH+ concerning the ordering o 



le ^S and '^A states. Therefore, we use 



MOLPRO along with state-of-the-art basis sets [38[ to compute several electronic states of 
these molecules in addition to the ^A states of interest. 

Diatomic molecules belong to the C200 point group; however, MOLPRO only uses Abelian 
point groups, necessitating the use of the C2V point group for our calculations. In this point 
group, the projection of electronic angular momentum along the diatomic axis (A) will 
transform as the following: A = 2n will transform as Al symmetry, A = — 2n as A2, 
A = 2n -|- 1 as Bl and A = —{2n + 1) as B2. As a result, the full A symmetry transforms as 
A1©A2, requiring the MOLPRO calculations to include both the Al and A2 symmetries. 



26 



27l | calculation at 



The MRCI ^ \m, calculations are performed after a MCSCF 
various internuclear separations (lA < r < lA). The term "active space" describes the set 
of orbitals being used to construct the configurations, as in Eq. p8|l . The active space is 
increased in one of two ways. One is to allow the electrons to reside in higher-lying orbitals 
that are nominally unoccupied. The second is to promote electrons from previously closed 
shells into these same orbitals. Both of these options allow more parameters with which to 
optimize the total energy, but it comes at the expense of an increased computational burden. 

11 



The basis sets and active space sizes are bench-marked against the dissociation energies 
of lli 18, bj. Our final choice was the ECP60MWB basis set of the Stuttgart/Cologne 



group 



32l | for the heavy atom along with t 



le aug-cc-pVTZ basis [39] for the hydrogen. 
38l | group stand for neutral atom quasi- 



m, 



The MWB bases of the Stuttgart/Cologne 
relativistic potentials. To make calculations manageable, these basis sets deal explicitly 
only with electrons outside of a 60-electron "core." The effect of these core electrons is 
incorporated through an effective core potential (ECP), which is also provided along with 
the basis sets. The ECP's include angular-averaged (i.e., scalar) relativistic effects. 

The active space chosen is of the form (Al,Bl,B2,A2) = (6,2,2,l) for both PtH+ and HfH+. 
This notation means that the first 6 Al, 2 Bl, 2 B2 and 1 A2 orbitals above the closed space 
are active. In PtH+ and HfH+, we close off orbitals in (2,1,1,0). Closed orbitals participate 
in optimizing the energy subject to the constraint that they are always doubly occupied. 
Thus, they do not participate in correlations. These calculations resulted in dissociation 
energies of 3.14eV(PtH^) and 2.49eV(HfH^) compared to previous theoretical energies of 
2.87eV and 2.50eV |l2i]. Dyall [1^ obtains a dissociation energy for PtH+ of 2.00eV while 
Zurita et al. jl9| obtain 2.22eV, both of which are significantly smaller than our result. Our 
deeper, and presumably more accurate, dissociation energies are most likely due to increased 
computational efficiency over the past decade. Calculations on PtH+ were completed in 
approximately 6 cpu-hours per point on a 2.4 GHz processor, while for HfH"*" they were 
completed in approximately 20 cpu- minutes per point. 

In order to further test the basis sets, we performed calculations on Pt and Hf to find 
the ionization energies. We ran a MCSCF/MRCI calculation on Pt, Pt^ and Pt^+ as well 
as for Hf, Hf"*" and Hf^+. We found the first and second ionization energies of Pt to be 
65000cm~^ and 143000cm~^ and of Hf to be 49600cm~^ and 115500cm"^. They agree with 



the empirical values |43l. l44| to within 6000 cm~^. 

Born-Oppenheimer potential curves for several symmetries are presented in Fig.(p]). The 
curves are generated by fitting an extended Rydberg function to the points obtained from 
the MOLPRO calculation, using the method of Ref. j^. All curves shown dissociate to 
the ground ^S state of H and the ^D state of the transition metal ion. To assign a Hund's 
case (a) label to a curve, we consult the MOLPRO output and assign a value based on the 
CI eigenvector. In turn, the eigenvector tells in which orbitals the electrons lie. There are 
several different configurations present, but the dominant configuration has a coefficient ho 

12 
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FIG. 1: Potential energy curves for PtH+ (a) and HfH+ (b), as obtained from MOLPRO output. 

very close to unity. As in the case of a FtH"*", we find there are several configurations, each 
having two unpaired electrons in the Al symmetry group. Upon examining the orbitals 
where these electrons lie, we find one has angular momentum projection A = 0, while the 
other has A = +2. Hence, this is a {a5) configuration. Some CI vectors have the same 
unpaired electrons but have paired electrons in a different orbital and/or symmetry group. 
In our approximation, these paired electrons play no role, and therefore they contribute 
nothing to the effective electric field within the molecule. Lastly, if this configuration is 
truly a '^A, there should be an equivalent energy configuration with one orbital in symmetry 
Al and another in symmetry A2. This is indeed what we find. 

Recall that Hund's case (a) describes a diatomic molecule with spin and orbital motion 
strongly coupled to the molecular axis. Hence the electronic part of the Hamiltonian is 
dominant followed by spin-orbit and then rotational coupling. Within the case (a) MRCI 
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calculation, both molecules nominally possess ^E ground states. However, both have low- 
lying ^A states. The energy difference Esj^ — Ei-^ is 1500 cm~^ for PtH+ and 2100 cm~^ 
for HfH"^. At this level we therefore agree qualitatively with the ordering of states given 

fH+. The equilibrium separations compare favorably to 



by Ref. |17| for PtH+, but not 



previous results for PtH"*" 



Q 



la 



19| and HfH"*" |l7|. Rotational constants were not quoted 



in these papers. In the ^A states of PtH+ and HfH+ we find r = 1.51A, Bg = 7.37cm ^ and 
r = 1.82Aand B = S.llcm'^ 

C. Estimating diagonal spin-orbit contributions 

Roughly speaking, the molecules inherit a large spin-orbit interaction from the heavy ion. 
Moreover, the spin-orbit constant A in PtH+ is expected to be negative, as it is in the Pt^ 
ion, thus shifting the ^A^ state lower in energy. For this reason it was previously anticipated 
that ^Aa is likely to be the absolute ground state of PtH+ [18]. By contrast, in HfH"*", A 
is expected to be positive and hence the ^Ai state is shifted down in energy. To estimate 
the size of the shift, we require knowledge of A in both molecules. Using methods outlined 
in Ref. [20], we estimate the value of A using van Vleck's "pure precession hypothesis" 
(42| . Under this approximation, the spin-orbit coupling arises from the heavy atom, and the 
atomic orbital angular momentum / of each electron is taken to be a good quantum number. 

Using this assumption, the "stretched state" ^Aa is represented by the following Slater 
determinant: 

I^Ag) = |ci5"sa"|. (23) 

In this expression, we make explicit the orbital angular momentum quantum numbers d 
and s] the index a stands for an electron spin aligned along the molecular axis, a = -1-1/2. 
The other states can be generated from this one by the application of suitable spin-lowering 
operators. We can now act upon wave functions of this type with the single-electron-spin- 
orbit operator, by recasting it in terms of operators acting on the individual electrons i: 

H^^ = AL-S 

= ^a4i-Si. (24) 

i 

Details of this procedure are presented in the Appendix. A more thorough treatment of 
spin-orbit and rotational coupling will be performed elsewhere J45[. 
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FIG. 2: Hund's case (c) potential curves for PtH^ (a) and HfH^ (b), labeled by the quantum 
number J7. These curves are obtained from those in Fig.lpQ) by adding spin-orbit corrections 
perturbatively, as described in the text. 



The effect of the spin-orbit mixing is to break the degeneracy of the different Q values 
and produce a revised set of potential curves, indexed by Hund's case (c) quantum numbers. 
These curves are shown in Fig.Q. Performing these calculations yields values oi A = 
— 1680cm~^ and A = 610cm~^ for the molecular spin-orbit constants of PtH"*" and HfH"*", 
respectively. This shift is more than enough to push the ^As curve below that of the ^S in 
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PtH"*" (Fig.(j2ti'))- Recall that the ^E state shifts little in energy because there is no spin or 
orbital motion, and therefore any shift in energy must be a second-order effect. The lowest 
fi = 3 curve is solely ^A^. Additionally, the ^A^ state of PtH+ is well described in Hund's 
case (a); there are no other states with fi = 3 nearby in energy (the closest we estimate is 
a ^$3 state ^ 30,000 cm^^ away that dissociates to the ^F limit of Pt"*'). 

The same results applied to HfH'*' tell a different story (Fig.(j2l3)). The value A is not 
sufficiently large to shift the ^Ai state below the ^Sq state. The energy difference between 
Q = 1 and i7 = states is 600 cm~^. There are other relativistic effects, such as the second- 
order correction to the spin-orbit constant and spin-spin and spin-rotation interactions, 
that can further push down the ^Ai state. For example, we can estimate the second- 
order-spin-spin parameter in terms of the atomic spin-orbit constant in the pure-precession 
approximation J30|. In these heavy elements it is the more dominant contribution to spin- 
spin. We find 

|(3A2|//.o|iA2p 



A(2) 



E:i/\ — El 

1 a2 



2 Es^ — Ei/\ 
^ 60cm"^ (25) 

These contribute to an uncertainty of the an order of magnitude smaller than the energy 
difference in the spin-orbit-corrected curves. It is conceivable that the sum total of these 
other relativistic effects, acting more strongly on the '^A state than on the ^S state, will 
reverse the ordering in Fig.(j2l3)- In principle, there is an uncertainty in the computed 
electronic energy using MOLPRO, but this is not trivial to estimate. States more subject 
to relativistic effects (like '^A) are less accurately computed compared to a ^S state. More 
detailed calculations — or, even better, experiments — will be required to decide this issue. 

If it turns out that the ground state of HfH+ is ^S, then the ^Ai state is metastable, 
and an estimate of its lifetime is desirable. On symmetry grounds we see that this lifetime 
should be very long, because of the fact that the transition would require AS* = 1, AQ = 1 
and AA = 2. J is the total angular momentum about the internuclear axis, J = L + S + N, 
where N is the mechanical rotation of the nuclei about their center of mass. The S and 
Q changes are are not dipole allowed in the case (a) basis. However, rotational couplings 
do not preserve the value of Q, and spin-orbit couplings do not preserve the value of S. 
Because of such interactions, the molecule is not a pure Hund's case (a) molecule. The 
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^A states may therefore by contaminated by small amplitudes of other case (a) states that 
are dipole allowed. These are very small effects and is the reason that these terms are 
usually neglected in the Born-Oppenheimer approximation. For example, we can compute 
the amount of contamination the ^IIi in the ^Ai state via spin-orbit. This was done to 
produce the curves in Fig.Q. Using the eigenvectors of the case (c) fi = 1 functions we 
compute the lifetime to be {\n = 1) = r/j^Ai) + (l^Ui) + ^I^Hi)): 

r ^ ^ea'{Es^, - Ei^^f 
~ 10 s. 

We see that lifetime is on the order of seconds. 

D. Estimating e^ and e^; the effective electric field 

As we have seen, the electric field inside a polar diatomic molecule is independent of 
the applied electric field, provided that the molecule is polarized. There is a relativistic 
enhancement factor determined from properties of the heavy atom and a projection of the 
electron spin, a, onto the molecular axis. Now we must determine, according to our working 
formula Eq. ()22|) . how much the s and p atomic orbitals are mixed because of the presence 
of the other atom. 

The values of es and ep are obtained from the molecular structure calculations, as de- 
scribed above. Wave functions generated for the a molecular orbital (not to be confused 
with the electron-spin projection a) and presented in Fig.Q- In each panel, the electronic 
wave function on the molecular axis is plotted as a function of z, the electron's coordinate 
along this axis (heavy line). For comparison, the s, p, and hydrogen atomic orbitals are 
also shown below the molecular orbital wave function. All wave functions are normalized 
to unity. Notice that the HfH"*" molecular ion in Fig.(jnb) is not as strongly mixed as in the 
case of PtH"*". This small mixing manifests itself as a smaller effective electric field for HfH"*" 
(not just due to the Z difference of 78 to 72). Qualitatively, this asymmetry illustrates the 
degree of s-p mixing that is vital to the observed electron EDM shift. By taking appropriate 
overlaps of atomic and molecular wave functions, we construct the effective electric field as 
defined in Sec. IIII Al 

To check the validity of this procedure we have performed the calculations of the expected 
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FIG. 3: Electronic wave functions (solid lines) for the relevant a molecular orbitals in the PtH^ 
(a) and HfH"^ (b) ions. In each figure, the lower panel shows the dominant s and p orbitals on the 
heavy ion, plus the s orbital on the hydrogen atom. All wave functions are normalized to unity. By 
projecting the molecular orbital onto the atomic orbitals, we find e^ = 0.75, Cp = 0.19 for PtII+, 
and es = 0.79, Cp = 0.09 for HfH+. 

electric field in molecules for which F^s has already been published; e.g., the ^S fluorides 
as well as the group VII hydrides and PbO in its metastable a^S state. Values of R 
(the internuclear separation) and term symbols are taken from the previous work on these 



molecules. Please see Ref. 
Ref. y for HrBr+, Refs. 

sets of the Stuttgart/Colo. 



Ref. m for Br, Ref. 



• m 




Q 



for BaF, Ref. |10| for PbF and HgF, Refs. |lj, [ij for HI+, 
for YbF, and Ref. [q] for PbO. We once again used the basis 



group for the heavy e 



for I, Ref. 



m for Yb, Ref. 



ements. Please see Ref. j35|] for Ba, 
for Pb and Ref. m for Pt, Hf 



and Hg. The lighter elements H and F use the aug-cc-pVTZ basis [39]. These results are 
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TABLE I: Comparisons of published values of -Fefr to our results using nonrelativistic software 



Molecule 


Published Fes (GV/cm) 


This work Feff (GV/cm) 


BaF 


7.r 


5.1 


YbF 


26^ 


43 


HgF 


99 -= 


68 


PbF 


-29'= 


-36.6 


PbO 


6.1'^ 


3.2 


HBr+ 


-0.02^ 


0.16 


HI+ 


-0.1^0.34/ 


0.57 


PtH+ 


N/A 


73 


HfH+ 


N/A 


-17 


HfF+ 


N/A 


-18 



'^Ref. 
''Ref. 
'^Ref. 
-^Ref. 
'^Ref. 
^Ref. 



i2\, "Ionic". 



presented in Table HI Although the literature values vary by nearly two orders of magnitude, 
our present method tracks most of them to within a factor of two. 

The lone exception is HBr^. The sign discrepancy between our value and that in Ref. 
(l2J | is analogous to the sign difference in HI+, as pointed out in Ref. [l^. Differences in 
the magnitude are likely due to Ref. (l2]'s use of an "ionic" bond, whereas HBr"*" is neither 
purely ionic nor covalent. Our method assumes nothing about the nature of the bond; 
rather, it allows the molecular structure software to calculate the electron distribution in 
the molecule and then projects the distribution onto atomic orbitals of the heavy atom. 

The values we obtain for the effective electric fields are: 73 GV/cm in the ^A^ state of 
PtH+ and -17 GV/cm in the ^Ai state of HfH+. Strikingly, these estimates are competitive 
with the largest ones previously identified. This result points to the imperative need for a 
complete electronic structure calculation to refine the predictions for these ions. 

Finally, at the suggestion of E. Cornell's group, we have used the same methods to 
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estimate the effective electric field for HfF"*". Like HfH"*", its ground state is a toss-up 
between ^Ai and ^Sq, and its effective electric field is F^s ~ —18 GV/cm. 

E. r^-doubling and critical field 

There remains the question of how easily polarized are the molecular states we have 
identified. The energy difference between states of differing parity Eq. (jl3|) in a A state is a 
fourth order effect in perturbation theory provided the '^^'^^Yl and ^•S'+i j] electronic states are 



fairly well separated in energy from the ^■S'+i^ state. Brown, et al. 4j| have used this idea 
to give formulas for estimating the effect in a manner consistent with the ideas of MuUiken 
and Christy 40], who worked out the effect in ^11 states. Thus, the fi-doubling is relegated 
to finding the A-doubling (the difference in energy between parity states comprised of A = 2 
and A = -2). 

In a ^A state, there are three parameters contributing to the energy splitting designated 
oa, Pa, and q/\. The parameters have the following forms 

5a = ^^^ ^ftY" { ) T^T—^. (26) 



AE1AE2AE: 




PA = 4V35— — — -— (27) 



^A - 8-— —-———, (28j 



AE1AE2AE3 
'VbYbYb 

'AE1AE2AE3 

where Va{b) is the matrix element of the spin-orbit (rotational) Hamiltonian between differ- 
ing states. (^Ai|AL_|'^ni) is an example of Va, whereas (^Ai|i?L_|^no) is an example of 
Vb- The AEi are the energy differences between the state of interest and the intermediate 
states. It is now straight-forward to estimate the order of magnitude of Edoub in the ^A 
states. 

In HfII+, we are concerned with the fi = 1 state, governed by the oa parameter. We find 
Edoub to be (recall that J = 1 in the ground state) 

ElZl = 25aJ(J+1) 

'^'' J(J + 1) 



(-EsA — -E3n)^(-E3A — Eij^ 
5.1 X 10"^ cm-^ 
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This Edouh leads to an Fcrit ~ IV/cm using the value of dm = 2.92 Debye obtained via the 
MOLPRO calculations. PtH+ is slightly more complicated. The parameter governing its 
splitting is q^ and connects the ^ Ai state to the ^ A3 state. So, we first find the value of the 
parameter and then use second-order perturbation theory within the ^A^ manifold. Thus 
we find (with J = 3 in the ground state) that 

j,n=3 ^ 2g^o( J( J + f ))3( J( J + 1) - 2)( J( J + 1) - 6) 

doub ^4^ _ g^^2 

4 ■ 10^51^2 



{4A - 85)2(E3A - ^3n)6(E3A - ^is)3 
3.32 X 10-2° cm"^ 



Using dm = 0.77 Debye in PtH+, from the MOLPRO calculation, we find Fcru ~ W^^V/cm, 
which is negligible for any purpose. 

IV. CONCLUSIONS 

The ^ A states of PtH+ and HfH+ satisfy our two most important criteria for experimental 
searches of the electron EDM. Both are easily polarized and expected to generate very high 
effective internal electric fields. However, several other experimentally relevant concerns 
remain to be resolved for these molecules, including: (1) whether ^Ai is the true ground 
state of HfH+; (2) at what laser wavelengths the molecules can be probed and detected; 
and (3) how their magnetic properties, e.g., g- factors, are influenced by channel couplings 
engendered by electric and magnetic fields and terms neglected in the Born-Oppenheimer 
approximation; i.e., intermediate coupling via J'^S", etc. These items demand further 
detailed study, and will be the subject of future work. 

We have endeavored to articulate a strategy for estimating the influence of the electron 
EDM on molecular structure, using perturbation theory and nonrelativistic electronic struc- 
ture calculation. This strategy has proved quite successful and should be useful to anyone 
wishing to assess the viability of a candidate molecule in the future. This estimate does not, 
of course, take the place of a complete relativistic many-electron calculation, which must be 
done as a follow-up. 
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APPENDIX A: SPIN-ORBIT IN VAN VLECK APPROXIMATION 

First we calculate in terms of Hund's case (a) functions and then in terms of the single 
electron operators. We present the wave functions in terms of Hund's case (a) numbers. 
These are presented by |A; 5*, S; J, fi) with J > Q. In terms of Slater determinants, the 
molecular wave functions are: 

\^A,) = |2;1,1;J,3) (Al) 

l^'Aa) = |2;1,0;J,2) (A2) 

= ^(|rf5"sa'^| + |d(5^sa"|) 
v2 

|=^Ai) = |2;1,-1;J,1) (A3) 



IdSf^sa^ 



1%) = |1;1,1;J,2) (A4) 

= |a7r S(J \. 

Once again normalizations are assumed in the determinants. The spin-orbit energy shift for 
the ^As is then 

(^AalAL.^.I^Ag) = 2A (A5) 

{\d5''sa''\\a{llsl + llsl)\\d5''sa'') = a(2(l/2) + 0(1/2)) 

= a. (A6) 

Now we apply the total spin lowering operator (given by S- = '^- s~) to Eq. ljAljl to obtain 
the ^A2 state. Since the projection of the spin, E, is zero, there is no shift. Therefore, the 
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molecular spin orbital constant A as half the atomic spin-orbit constant a. The spin-orbit 
parameter is the same for all states within the ^A manifold. 

There are higher-order effects via coupling to other fi = 3 electronic states. Yet, the 
only other i7 = 3 states come from ^$ or ^$ states, which lie far off in energy. They can 
be evaluated to see the second-order contribution to the molecular spin-orbit constant via a 
calculation similar to the one below describing the coupling oi Q = 2 states. 

Off-diagonal effects of spin orbit can be evaluated using the same one-electron operators. 
First we write the off diagonal elements using spherical tensor notation. 

^ E ^''i' + ^2S2'- (AT) 

Then we evaluate the coupling between Eq. ()A2|l and Eq. ()A4|) . The only one electron element 
that works is l^s^ to take ^112 to ^A2. Only including the terms that contribute, we find 

{J^^lj'^rll*-.^"!) = ^V(6-2)(3/4+l/4) 

- ^2 <^«' 

These calculations take into account the permutations of the indices within the deter- 
minant. This is what cancels the normalization. We have ignored the effects of the other 
electrons in the Slater determinant that belong to paired orbitals. These can, in principle, 
contribute via permutations of the indices and raising/lowering the angular momentum. 
These types of contributions will be evaluated in future studies. 
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